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摘 要 : 部 善 库 木 塔 格 沙漠 是 一 个 与 城市 相连 的 沙漠 ,蕴含 着 丰富 独特 的 微生物 菌株 资源 。 为 了 探究 部 善 库 木 塔 
格 沙 漠 的 微生物 群落 结构 以 及 其 与 环境 因子 之 间 的 关系 ,采集 了 部 善 库 木 塔 格 沙漠 北 缘 4 个 样 点 的 表层 (0~5 em) 
和 深层 (45~50 cm) 的 沙土 样本 ,并 测定 了 其 含水 量 (WC) .pH 总 碳 (TC) 、 总 氮 (TN) 总 有 机 碳 (TOC) 和 电导 率 
CEC) ;利用 高 通 量 测 序 技术 ,研究 不 同 采样 点 .不 同 深度 沙土 样本 的 细菌 群落 结构 ;运用 Spearman 相关 性 分 析 环 境 
因子 和 细菌 群落 结构 之 间 的 相关 性 ,探讨 影响 细菌 群落 结构 的 主要 因素 。 结 果 表 明 :部 善 库 木 塔 格 沙 漠 北 缘 沙 十 


样本 呈 碱 性 ; 门 水 平 上 的 优势 菌 是 放 线 菌 (Actinobacteria) 、 厚 壁 菌 (Firmicutes) 和 变形 菌 (Proteobacteria) , 属 水 平 上 


的 优势 菌 是 unclassified_f_Micrococcaceae , Bacillus , Sphingomonas ,Escherichia-Shigella, PCoA 分 析 表 明 ,部 善 库 木 


塔 格 沙漠 北 缘 不 同样 点 之 间 的 细菌 群落 结构 存在 显著 差异 (P < 0.05) ;表层 和 深层 样本 的 物种 多 样 性 没有 显著 差 
异 。RDA 分 析 表 明 ,WC、TC、TN 和 pH 显著 影响 表层 样本 的 细菌 群落 结构 (P<0.05) ,TOC 含 量 显 著 影 响 次 层 样本 
的 细菌 群落 结构 (P< 0.05)。 本 研究 为 后 期 部 善 库 木 塔 格 沙漠 及 相似 生境 中 微生物 资源 的 发 掘 提供 了 理论 依据 。 
关键 词 : 部 善 库 木 塔 格 沙漠 ; 高 通 量 测序 ， 细菌 群落 结构 ; 环境 因子 ; 理化 性 质 


由 于 沙漠 高 辐射 .昼夜 温差 大 、 寡 营养 等 特点 ， 
使 得 沙漠 中 能 生存 的 动 植 物 极 少 ,因此 微生物 成 为 
了 沙漠 中 的 重要 组 成 部 分 ”"。 痢 善 库 木 塔 格 沙漠 位 
于 新 疆 东 部 ,部 善 县 以 南 (91"25~94*00' 眉 ,38"20~ 
40°50'N) ,面积 约 为 24243 km?, tH T HR Es Ep o EL 
JUR 2 km, 被 称 为 “城中 沙漠 "中, 具有 “ 绿 不 退 , 沙 不 
进 ” 的 特殊 地 理 环境 ,其 特殊 的 生境 缠 含 着 独特 的 
微生物 群落 。 目 前 ,对 于 部 善 库 木 塔 格 沙漠 微生物 
方面 的 研究 很 少 ,少见 的 报道 是 关于 从 该 沙漠 中 分 
离 到 的 放 线 菌 菌株 产 淀粉 酶 的 研究 ”。 

传统 的 微生物 培养 技术 由 于 培养 基 成 分 .接种 
量 .培养 温度 .培养 时 间 等 的 限制 ,可 培养 出 的 微 生 
物 总 量 不 足 10% , 极 大 的 限制 了 人 们 对 整个 微生物 
群落 结构 的 了 解 "。 高 通 量 测序 由 于 准确 性 高 妆 
据 量 大 、 价 格 成 本 低 等 特点 ,在 当今 科研 领域 被 广 
泛 使 用 5” 。 基 于 高 通 量 测序 技术 的 研究 表明 , 沙漠 
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中 的 细菌 优势 菌 门 主要 是 放 线 菌 门 . 厚 壁 菌 门 AE 
PAY . 拟 杆菌 门 酸 杆菌 门 AA PRECES ERIT], 
然而 ,目前 对 于 部 善 库 木 塔 格 沙漠 微生物 群落 结构 
的 研究 鲜 见 报道 。 

土壤 特性 和 环境 因子 对 土壤 微生物 群落 的 结 
FA RACE Me), VD GEE AS RSE Sit FES FRY RF 
性 使 得 沙子 土壤 的 理化 性 质 对 沙漠 微生物 群落 的 
结构 影响 较 大 ,有 研究 表明 ,盐分 是 影响 古 尔 班 通 
古 特 沙漠 微生物 群落 结构 的 决定 因素 ”, 而 总 磷 、 
SO; 和 HCO; 含量 对 塔克拉玛干 沙漠 东 缘 的 细菌 群 
落 结构 具有 显著 的 影响 呈 。 研 究 环境 因子 对 沙漠 
微生物 群落 结构 的 驱动 ,对 了 解 沙 漠 生 态 系统 的 脆 
弱 性 至 关 重 要 。 

本 研究 利用 16S rRNA 基因 高 通 量 测序 技术 , 通 
过 对 部 善 库 木 塔 格 沙漠 北 缘 沙 土 样本 进行 细菌 群 
落 结 构 以 及 环境 因子 对 细菌 群落 结构 的 影响 进行 
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研究 ,以 期 为 后 期 部 善 库 木 塔 格 沙漠 的 微生物 研究 
提供 理论 依据 。 


1 材料 与 方法 


1.1 样品 采集 

样 点 选取 于 部 善 库 木 塔 格 沙漠 北 缘 无 人 为 干 
扰 和 无 植被 的 区 域 (图 1) ,样本 采集 于 2018 年 9 月 ， 
共 4 个 采样 点 (Kl1、K2、K3、K4), 每 个 采样 点 分 别 采 
集 表 层 (S:0~5 cm) 和 深层 (D:45~50 em) HA T 
度 。 表 层 样本 佩戴 无 菌 手套 直接 采集 ,次 层 样本 先 
使 用 铁 鳞 控 取 并 用 卷 尺 测量 45~50 cm 深度 后 ,佩戴 
无 菌 手套 采集 样本 。 每 个 采样 点 按照 五 点 法 分 别 
采集 3 个 平行 样本 ,3 个 平行 样本 间距 约 50 cm ,每 
个 样本 采集 约 100 g, 共 采集 24 份 样本 。 采 集 的 样 
本 装 于 无 菌 密 封 采样 袋 中 ,低温 运输 于 实验 室 ， 
于 -20 CHEATER. 
1.2 样品 理化 性 质 的 测定 

样品 含水 量 (WC) 通 过 烘 干 法 进行 测定 '”。pH 
样品 按照 去 离子 水 为 1:2.5(W/V) 混 匀 静 置 30 min 
后 用 pH 计 进 行 测定 中。 电导 率 (EC) 样 品 按照 去 离 
子 水 为 1:5 (W/V) 混 匀 静 置 30 min 后 用 电导 仪 
(DDSJ-308A, 上 海 雷 磁 ) 进 行 测定 后 。 样 品 的 总 碳 
(TC) 总 有 机 碳 (TOC) 和 总 氮 (TN) 使 用 元 素 分 析 仪 
(Elementar Vario-EL , (£ Eg] ) JEfT jill qz ^, 
1.3 样品 总 DNA 的 提取 和 16S rRNA 基因 扩 增 

称 取 0.5 g 样 品 ,使 用 土壤 试剂 盒 (Qiagen , 德 
) 对 样品 的 总 DNA 进行 提取 。 吸 取 1 kL 提取 好 
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的 样本 DNA, 使 用 DNA 浓度 测定 仪 (NanoDrop 
2000) 测 定 其 浓度 。 选 取 16S rRNA 基因 的 V3~V4 
区 进行 扩 增 ,所 用 的 特异 性 引物 是 338F-806R"*。 
PCR 扩 增 条 件 为 :94 "C UAE TE 5 min;95 CC 变性 30 
s;55 % 退 火 30 s;72 CHE f 45 s;72 % 后 延伸 10 
min, 共 30 个 循环 。 将 PCR 扩 增产 物 浓 度 大 于 0.5 
ng* ML 的 样本 进行 扩 增 子 测序 。 
1.4 扩 增 子 测序 及 数据 分 析 

样本 的 16S rRNA V3~V4 区 扩 增 产物 在 Tllumi- 
na Miseq PE300 平台 上 进行 测定 。 通 过 QIME 软件 
包 对 原始 数据 进行 质量 过 滤 ' ;使 用 软件 fastp 和 
FLASH 对 序列 数据 进行 去 杂 ; 去 杂 后 的 序列 以 97% 
的 相似 度 聚 类 为 操作 分 类 单元 (OTUs) ;使 用 Uparse 
软件 对 OTU 进行 聚 类 ,运用 了 USEARCH11-uparse 
算法 ;采用 RDP classifier 软件 对 97% 相 似 水 平 的 
OTU 代表 序列 进行 分 类 学 比 对 ;物种 由 SILVA138 
细菌 数据 库 进 行 注 释 。 

采用 SPSS 软件 (version 23.0) 对 沙子 土壤 样本 
的 环境 因子 和 Alpha 多 样 性 指数 进行 统计 分 析 ,使 
用 单 因 素 方 差分 析 进 行 差异 显著 性 检验 。 利 用 
Mothur 计算 不 同 随机 抽样 下 的 Alpha 多 样 性 指数 ， 
再 使 用 软件 Origin 2021 绘制 稀释 曲线 图 。 运 用 R 语 
言 (版 本 3.3.1) 的 vegan 包 进行 元 余 分 析 (RDA) ,用 
来 反映 不 同样 点 的 环境 因子 与 细菌 群落 之 间 的 相 
关 性 。 使 用 RR 语言 ( 版 本 3.3.1) 的 heatmap 包 绘制 
heatmap 图 。 利 用 Spearman 相关 系数 显示 不 同样 点 
的 环境 因子 与 不 同门 之 间 的 相关 性 。 基 于 Bray- 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2020)4634 号 ,对 底 图 边界 无 修改 。 


1 部 善 库 木 塔 格 沙漠 北 缘 采 样 点 示意 图 


Fig. 1 Sampling map of the northern margin of Shanshan Kumtag Desert 
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Curtis 距离 的 主 坐 标 分析 (PCoA ) 用 来 分 析 不 同样 点 
间 的 细菌 群落 相似 性 和 差异 性 。LEfSe 软件 用 来 分 
析 所 有 表层 和 深层 样 点 之 间 的 差异 物种 ,LDA E 
设置 为 2。 


2 结果 与 分 析 


样本 的 理化 性 质 
沙子 土壤 样本 环境 因子 分 析 表 明 ,各 样 点 之 间 
的 理化 性 质 差 异 较 大 ,总 体 来 说 沙土 样本 呈 碱 性 
( 表 1)。 表 层 不 同样 本 的 pH EC TC TOC 和 WC 均 
星 显 著 差异 (P<0.05), 只 有 TN 不 显著 ;深层 不 同样 
本 除了 EC 不 显著 ,pH、TC、TN .TOC 和 WC 均 呈 显著 
差异 。 
2.2 细菌 群落 组 成 

对 部 善 库 木 塔 格 沙漠 北 缘 24 个 沙土 样本 测序 
后 进行 分 析 , 双 端 测序 共 获 得 2585300 个 序列 ,获得 


2.1 


的 总 碱 基 数目 为 1556350600。 双 端 序 列 拼接 质 控 
优化 后 的 序列 是 1292650 条 ,优化 后 的 碱 基数 目 是 
563481141, OTU 在 97% 相 似 度 水 平 下 聚 类 后 得 到 
2450 个 OTU ,在 SILVA138 数 据 库 中 注释 后 共 得 到 
48 个 门 ,125 个 纲 ,287 个 日 ,466 个 科 ,867 个 属 ， 
1378 个 种 。 

部 善 库 木 塔 格 沙 江北 缘 不 同样 点 的 微生物 群 
落 组 成 存在 一 定 差 异 ,将 平均 相对 丰 度 小 于 1% 的 
类 群 归 为 其 他 (others)。 按 照相 对 丰 度 的 高 低 依次 
排序 ,在 门 水 平 相对 丰 度 大 于 1% 的 菌 门 是 放 线 菌 
门 (Actinobacteria) 、 变 形 菌 门 (Proteobacteria) , EE fal 
门 (Firmicutes)、 拟 杆菌 门 (Bacteroidetes)、 芽 单 胞 菌 门 
(Gemmatimonadetes), Z& ?5 [al | ] (Chloroflexi), Deino- 
coccota、Patescibacteria( 图 2)。 整 体 来 说 ,无 论 是 表 
层 样 本 还 是 深层 样本 ,优势 菌 门 (相对 丰 度 > 10% ) 
是 一 样 的 ,都 是 放 线 菌 门 变形 阔 门 和 厚 壁 菌 门 。 
K2 和 K3 样 点 的 放 线 菌 门 相对 丰 度 (47.3%~61.1%) 


表 1 不 同样 点 样本 的 理化 性 质 


Tab.1 Physical and chemical properties of samples at different points 


样本 名 称 WC pH EC/(mS - cm ?) TN/% TC/% TOC/% 
KI_I_S 0.64 9.46 150 0.040 0.670 0.281 
K125 1.30 9.49 143 0.035 0.745 0.321 
K135 1.55 9.67 171 0.044 0.732 0.300 
K218S 0.52 8.66 334 0.040 0.485 0.197 
K228 0.53 8.8 350 0.037 0.409 0.196 
K235 0.42 8.98 236 0.029 0.494 0.182 
K315 0.72 9.14 594 0.033 0.383 0.238 
K328 0.53 9.21 523 0.029 0.438 0.160 
K335 0.49 9.22 616 0.030 0.336 0.178 
K4 18 0.47 9.42 258 0.026 0.513 0.185 
K425 0.37 9.3 282 0.032 0.613 0.236 
K435 0.30 9.32 408 0.043 0.508 0.189 
KI 1D 0.28 8.89 663 0.034 0.616 0.235 
K1_2_D 0.35 8.83 633 0.037 0.516 0.231 
K1_3_D 0.26 8.97 278 0.035 0.688 0.247 
K21D 0.26 9.17 345 0.021 0.430 0.161 
K22D 0.35 9.06 332 0.027 0.456 0.214 
K23D 0.28 9.05 358 0.021 0.353 0.152 
K31D 0.34 9.29 406 0.029 0.527 0.139 
K32D 0.45 9.29 424 0.028 0.385 0.128 
K3_3_D 0.56 9.35 287 0.042 0.408 0.176 
K4 1D 0.25 8.76 334 0.034 0.571 0.215 
K42D 0.35 8.81 282 0.026 0.498 0.220 
K43D 0.37 8.64 466 0.037 0.510 0.211 


注 :EC 为 电导 率 ;TN WAR; TC 为 总 碳 ;TOC 为 总 有 机 碳 ; WC 为 含水 量 。 
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图 2 门 水 平 上 沙子 样本 细菌 群落 结构 


Fig. 2 Bacterial community structure of sand samples at phylum level 


明显 高 于 K1 A K4 BESS (26.8% ~44.9% ) ,而 K2 和 K3 
样 点 的 厚 壁 菌 门 的 相对 丰 度 (10.9%~13.4% ) 明 显 低 
T K1 和 K4 样 点 (20.9%~40.3% )。 变 形 菌 门 在 各 样 
点 的 相对 丰 度 相差 不 明显 。 尽 管 各 样 点 的 优势 菌 
门 相同 ,但 是 在 同一 个 样 点 中 ,表层 和 深层 的 物种 
相对 丰 度 是 有 差异 的 。 比 如 在 Kl1 样 点 深层 样本 中 
拟 杆菌 门 芽 单 胞 菌 门 和 Deinococcota 的 相对 丰 度 
分 别 是 2.8% .3.9% 和 2.4% ,而 在 K1 样 点 表层 中 的 相 
对 丰 度 分 别 是 0.7% .0.2% 和 0.4%。 绿 弯 菌 门 在 K4 
样 点 表层 中 的 相对 丰 度 为 1.3%% ,而 在 K4 样 点 深层 
样本 中 相对 丰 度 是 0.5%。 

在 属 水 平 ,相对 丰 度 ( > 1%) 依 次 从 高 到 低 分 别 
为 : unclassified_f__Micrococcaceae , Bacillus , Sphin- 
gomonas , Escherichia- Shigella , Microbacterium , Brevun- 


dimonas , Pseudomonas , Citrobacter , Acinetobacter , Sphin- 


1.0 


0.8 


属 水 平 相对 丰 度 


KIS KID K2s K2D K3S 
样品 采样 点 


gobium , Stenotrophomonas .norank { norank _o__Fran- 
kialesnorank, norank f  Longimicrobiaceae, Entero- 


coccus , Deinococcus , Romboutsia 和 Flavisolibacter (图 

)。 其 中 最 优势 的 属 为 unclassified_ f ^ Micrococcace- 
上 5%~53.9% ) „Bacillus (4.1%~38.5% ) ,Sphingomo- 
nas (5.5% ~21.3% ) 和 Escherichia-Shigella (17.4%). 7S 
同样 点 的 优势 属相 对 丰 度 存在 差异 。unclassi- 
fied_f_Micrococcaceae 的 相对 丰 度 表现 为 :K2> 
K3 > K4 > K1 ; Bacillus fll Escherichia-Shigella 的 相对 
E BE RMN : K1 > KA > K3>K2; Sphingomonas 的 
相对 丰 度 表现 为 :K3>K2>K4>Kl。norank f - 
norank o. Frankialesnorank, norank f  Longimicro- 
biaceae „Deinococcus 和 Flavisolibacter 在 Kl1 样 点 深层 
样本 中 的 相对 丰 度 (5.3% ) 明 显 高 于 其 余 样 点 ;Rom- 
boutsia 在 K2 样 点 表层 样本 中 的 相对 丰 度 最 高 。 同 


m unclassified f Micrococcaceae 
= Bacillus 
m Sphingomonas 
m Escherichia-Shigella 
m Microbacterium 
m Brevundimonas 
m Pseudomonas 
m Citrobacter 
m Acinetobacter 
| Sphingobium 
m Stenotrophomonas 
m norank f norank o Frankiales 
m norank f Longimicrobiaceae 
m Enterococcus 
m Deinococcus 
m Romboutsia 
m Flavisolibacter 
m others 


K3D K4S K4D 


图 3 属 水 平 上 沙子 样本 细菌 群落 结构 


Fig. 3 Bacterial community structure of sand samples at genus level 


202308.00666v1 


ChinaXiv 


1362 To 


ChinaXiv 合 作 期 刊 


区 研 R 40 卷 


一 样 点 中 ,表层 和 深层 在 属 水 平 上 物种 的 相对 丰 度 
也 存在 明显 差异 。Bacillus 在 Kl1 样 点 表层 样本 的 相 
对 丰 度 (38.5% ) 明 显 高 于 Kl 样 点 深层 样本 
(24.8% ) ; Microbacterium 在 K4 样 点 深层 样本 的 相对 
丰 度 (5.9%) 明 显 高 于 K4 样 点 表层 样本 (3.3%)。 
2.3 表层 和 深层 样本 物种 差异 分 析 

LEfSe 结果 表明 ,在 门 水 平 上 Campilobacterota 
和 NB1-j 在 表层 样本 中 显著 富 集 (P < 0.05) ,Sumerl- 
eota 在 深层 样本 中 显著 富 集 (P<0.05)。 属 水 平 上 
表层 显著 富 集 的 物种 有 17 个 (图 4) ,分别 是 Collin- 
Helicobacter , norank f norank o WCHBI- 41, 


norank f norank o Actinomarinales , norank, f noran- 


sella, 


k_o_Clostridia_vadinBB60_group , Prevotellaceae_UC- 


G-001 , Psychrobacter , norank_f_Bacteroidales_RF16_ 


group. Streptomonospora , Family_xlll_AD3011_group , 


Chelativorans , Parabacteroides , Prevotellaceae NK3B31 . 


group, Aerococcus , Fusobacterium , unclassified, f. Os- 
cillospiraceae 和 Pseudochrobactrum (P < 0.05) , 深层 


Pseudochrobactrum | 


unclassified f Oscillospiraceae : 
Fusobacterium : 
Aerococcus d 


显著 富 集 的 物种 有 14 个 ,分 别 是 norank_f_norank o 
Frankiales, norank f Longimicrobiaceae, 、 norank_f_ 
LWQ8, unclassified_f_Oxalobacteraceae , Rubellimicro- 
bium, Geodermatophilus , Noviherbaspirillum , norank_ 
f_norank_o_Saccharimonadales , unclassified f Nocar- 
dioidaceae , Oxalicibacterium , Ramlibacter , Ammoniphi- 
lus ,norank, f Symbiobacteraceae ^l Dietzia (P < 0.05 ) 
24 细菌 Alpha 多 样 性 和 Beta 多 样 性 

部 善 库 木 塔 格 沙漠 北 缘 细 菌 Alpha 多 样 性 指数 
统计 结果 ( 表 2) 显 示 , 样 品 的 覆盖 度 均 达到 99.9% 以 
上 ,证 明 测序 的 数据 可 靠 有 效 。 整 体 来 看 ,表层 和 
深层 之 间 的 样本 香农 指数 .辛普森 指数 Chao 1 指 
数 、 ACE 指数 和 系统 发 育 树 多 样 性 指数 没有 显著 差 
异 , 表 明 表 层 和 深层 样本 的 多 样 性 没有 明显 差异 。 
不 同样 点 之 间 的 Chao 1 指数 和 ACE 指数 存在 显著 
差异 (P<0.05), 而 香农 指数 、 辛 普 森 指数 和 系统 发 
育 树 多 样 性 指数 没有 显著 差异 ,表明 不 同样 点 之 间 
细菌 群落 的 丰富 度 明 显 不 同 ,而 多 样 性 没有 显著 差 


| 国 表层 


Prevotellaceae NK3B31 group [El 
Parabacteroides [SSS 
[eg —— 


Chelativorans 

Family xlll AD3011 group 
Streptomonospora 

norank f Bacteroidales RF16 group 
Psychrobacter 


Prevotellaceae UCG-001 [h] 
norank f norank o Clostridia vadinBB60 group [p] 


norank f norank o Actinomarinales 
norank f norank o WCHBI-41 
Helicobacter 

Collinsella 


norank f norank o Frankiales 
norank f Longimicrobiaceae 
norank f LWQ8 

unclassified f Oxalobacteraceae 
Rubellimicrobium 
Geodermatophilus 
Noviherbaspirillum 

norank f norank o Saccharimonadales 
unclassified f Nocardioidaceae 
Oxalicibacterium 

Ramlibacter 

Ammoniphilus 

norank f Symbiobacteraceae 
Dietzia 


| H : E r . - 
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10 05 00 05 10 
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LDA SCORE (lg) 


图 4 属 水 平 上 表层 和 深层 样本 的 物种 差异 


分 析 


Fig.4 Species difference analysis in Surface and deep Samples at genus level 
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表 2 不 同样 本 的 Alpha 多 样 性 指数 
Tab. 2 Alpha diversity index at different samples 


样本 名 称 香农 指数 辛普森 指数 ACE 指数 Chao 1 指数 敌 盖 度 /9% 系统 发 育 树 多 样 性 指数 
KL15 2.390 0.215 254.967 262.154 1.000 31.200 
K1 2 S 2.766 0.143 336.400 338.813 0.999 43.398 
K135 1.974 0.215 315.097 318.561 0.999 32.058 
K215 2.277 0.317 284.391 284.000 1.000 32.017 
K225 1.862 0.386 220.833 222.688 1.000 28.580 
K235 3.583 0.156 616.427 633.067 0.999 79.401 
K315 2.622 0.170 245.595 249.200 1.000 31.607 
K325 2.468 0.233 409.968 413.000 1.000 53.739 
K335 2.114 0.247 182.997 184.500 1.000 22.578 
K4 15 3.922 0.138 716.334 780.650 0.999 81.842 
K428 2.277 0.302 418.711 435.719 0.999 45.497 
K435 3.251 0.105 657.699 662.017 0.999 73.443 
KL 1D 2.856 0.258 504.457 508.081 0.999 51.449 
K1_2_D 3.986 0.068 638.214 643.000 0.999 57.332 
K1_3_D 2.827 0.155 579.787 599.192 0.998 48.243 
K2_1_D 1.843 0.425 230.223 230.273 1.000 31.763 
K22D 2.760 0.147 291.824 304.214 0.999 32.989 
K2 3_D 1.989 0.342 187.780 187.000 1.000 31.652 
K31D 2.350 0.345 389.315 407.083 1.000 59.362 
K32D 1.921 0.331 229.616 231.929 1.000 30.269 
K33 D 2.260 0.278 306.728 308.120 0.999 37.188 
K4 1 D 3.134 0.112 468.690 468.333 0.999 63.828 
K42 D 2.626 0.163 335.807 338.043 1.000 36.498 
K43 D 2.114 0.296 227.945 226.688 1.000 26.739 


e Kl 
a K2 
+ K3 
4 K4 


Fo PARE K4 ASE E EROS ,深层 样 点 K1 的 丰 
富 度 最 高 。 

利用 PCoA 方 法 ,分析 部 善 库 木 塔 格 沙漠 北 缘 
不 同样 点 之 间 细 菌 群落 的 相似 性 和 差异 性 (图 5)， 
PC1 轴 解 释 了 54.89% 的 差异 性 ,PC2 轴 解释 了 
13.4% 的 差异 性 ,整体 上 不 同样 点 的 细菌 群落 组 成 
存在 显著 差异 (P<0.05)。K1 和 K4 样 点 与 K2 和 K3 
样 点 相 比较 , 分 布 在 PC1 轴 的 样 点 离散 程度 高 ,说明 -04 -02 0 02 04 06 08 
K1 和 K4 样 点 的 组 内 样本 差异 较 大 。 
2.5 环境 因子 对 细菌 群落 结构 的 影响 图 5 不 同样 点 的 Beta 多 样 性 分 析 

为 了 分 析 环 境 因 子 对 部 善 库 木 塔 格 沙漠 北 缘 Fig. 5 Beta diversity analysis at different points 
AH bal AE Yes HS PA EGRE, ANE SESE RDA SPT TE CBR Yb IER ge RN Pe ETE AH CP < 0.05). OR 
OTU 水 平分 析 了 环境 因子 与 部 善 库 木 塔 格 沙漠 层 样本 中 (图 6b) ,第 一 轴 解 释 了 62.51% ,第 二 轴 解 
北 缘 细 菌 群落 之 间 的 相关 性 (图 6)。 结 果 表 明 , 采  ” 释 了 4.79%, 其 中 只 有 TOC 显著 影响 深层 样本 的 细 
样 深 度 对 环境 因子 与 细菌 群落 结构 的 影响 不 同 。 ARAETA CP < 0.05) 
表层 样本 中 (图 6a) ,第 一 轴 解 释 了 64.56% ,第 二 为 了 进一步 分 析 环 境 因子 对 不 同 物种 的 影响 ， 
MERET 9.41%, WC TC TN 和 p 卫 均 显 著 影 响 部 善 库 ”本 研究 选择 了 相对 丰 度 前 10 的 细菌 门 绘制 了 heat- 


PC2 (13.4%) 
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6 表层 (a) 和 深层 (b) 样 本 的 细菌 群落 与 环境 因子 之 间 的 相关 性 分 析 


Fig.6 Correlation analysis between bacterial communities and environmental factors in surface (a) and deep (b) samples 


(8) 相关 系数 
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(b) 相关 系数 
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图 7 表层 (a) 和 深层 (b) 样 本 在 门 水 平 上 的 不 同 物种 与 环境 因子 之 间 的 相关 性 分 析 


Fig. 7 Correlation analysis between different species and environmental factors in surface (a) and deep (b) samples at phylum level 


map 图 ,用 来 评估 物种 与 环境 变量 之 间 的 相关 性 (图 
7)。 结 果 表 明 :表层 样本 ,EC 5j Deinococcota 显著 负 
HÆ CP « 0.05) ;TC 与 Deinococcota 和 Firmicutes 显 
3E EAA (P « 0.05) ; TOC 5j Spirochaetota fib E fi 4H 
关 (P<0.05)。 对 于 深层 样本 ,WC 与 Deinococcota 
和 Chloroflexi 显著 负 相 关 (P<0.05) ; EC 与 Fir- 
micutes 显著 负 相 关 (P<0.05) ;TC 与 Bacteroidota、 
Deinococcota , Chloroflexi 和 Patescibacteria 显著 正 相 
X,+ Actinobacteriota 显著 负 相关 (P<0.05);TOC 与 
Gemmatimonadota , Deinococcota 和 Patescibacteria 显 


著 正 相关 ,与 Actinobacteriota 显著 负 相 关 (P< 0.05) 
3 讨论 


由 于 环境 的 异 质 性 ,驱动 不 同 深度 沙漠 微生物 
菌 群 变化 的 因素 也 就 不 同 。 本 研究 结果 表明 ,WC、 
TC、TN 和 pH 均 对 表层 样本 的 细菌 群落 结构 有 显著 
影响 ,而 深层 样本 的 细菌 群落 结构 仅 受 TOC 的 影 


响 。 由 于 表层 样本 直接 暴露 于 空气 中 , 受 外 界 极端 
环境 的 影响 较 大 , 而 微生物 生存 离 不 开水 分 、 碳 源 
和 氮 源 ,因此 ,WC、TC 和 TN 显著 影响 表层 微生物 菌 
群 的 结构 。 微 生物 利用 的 碳 源 主要 是 有 机 碳 , 因 此 
有 机 碳 的 含量 与 微生物 菌 群 的 代谢 活动 密切 相 
关 。 深 层 样本 的 采样 深度 是 45~50 om, 与 表层 相 比 
具有 较 稳 定 的 内 部 环境 ,外 界 环 境 变化 对 其 影响 较 
小 ,因此 ,深层 样本 中 微生物 菌 群 的 结构 主要 受 有 
机 碳 的 影响 。 痢 善 库 木 塔 格 沙漠 北 缘 不 同样 点 .不 
同 深度 的 样本 pH 均 旦 碱 性 ,这 与 塔克拉玛干 沙漠 
"和 古 尔 班 通 古 特 沙 漠 '" 的 研究 结果 相似 。 人 研究 
表明 ,pH 是 影响 土壤 细菌 群落 结构 的 重要 指标 。 
此 外 ,在 同一 样 点 不 同 深度 样本 的 土壤 理化 因子 也 
存在 差异 ,如 Ki1 样 点 表层 样本 的 WC 含量 最 高 ,而 
Kl 样 点 深层 样本 的 WC 并 不 是 最 高 ;EC 在 K3 样 点 
的 表层 样本 最 高 , 而 深层 不 同样 点 的 EC 含量 差异 
不 显著 。 
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本 研究 结果 表明 ,不同 环境 因子 对 表层 和 深层 
样本 的 优势 菌 门 的 影响 不 同 ,例如 WC 对 表层 样本 
的 酸 杆菌 门 影 响 显 著 , 而 对 深层 样本 的 酸 杆菌 门 影 
响 不 显著 ;pH 对 表层 样本 放 线 菌 门 影响 显著 ,而 对 
深层 样本 放 线 菌 门 影 响 不 显著 ;EC 对 深层 样本 厚 壁 
菌 门 影响 显著 ,而 对 表层 样本 厚 壁 菌 门 影响 不 显 
车 。 这 说 明 表 层 和 深层 样本 对 菌 群 结 构 具 有 自己 
独特 的 驱动 模式 。Bahadur 等 揭示 了 柴 达 木 沙 漠 
中 变形 菌 门 与 TOC 显著 相关 ,Zhang 等 ”对 腾 格 里 
沙漠 的 研究 具有 相同 结果 ,然而 本 研究 结果 表明 无 
论 是 表层 样本 还 是 深层 样本 ,变形 菌 门 与 TOC 相关 
性 不 显著 。 

Beta 多样 性 分 析 表 明 ,部 善 库 木 塔 格 沙漠 北 缘 
不 同样 点 的 细菌 群落 结构 存在 显著 差异 (P<0.05)， 
这 可 能 是 造成 不 同样 点 细菌 群落 丰富 度 差 异 的 原 
因 。 另 外 ,部 善 库 木 塔 格 沙漠 北 缘 表 层 和 深层 样本 
的 细菌 群落 多 样 性 和 丰富 度 都 没有 显著 差异 ,这 可 
能 是 由 于 采样 深度 不 馆 的 原因 。 不 同 生 境 沙 漠 定 
殖 的 微生物 优势 菌 门 存在 差异 ,人 研究 表明 沙漠 表皮 
存在 生物 结 皮 的 优势 菌 门 主要 是 蓝藻 门 ,比如 十 尔 
班 通 古 特 沙 漠 '"\ 腾 格 里 沙漠 和 莫 哈 维 沙漠 ”|。 
部 善 库 木 塔 格 沙漠 北 缘 无 论 是 表层 样本 还 是 深层 
样本 ,优势 菌 门 都 是 放 线 菌 门 .变形 菌 门 和 厚 壁 菌 
门 , 其 中 各 样 点 的 相对 丰 度 最 高 的 门 是 放 线 菌 门 ， 
这 与 塔克拉玛干 沙漠 "的 优势 菌 门 一 臻 。 放 线 菌 
门 能 产生 内 生 孢 子 用 来 抵御 极端 环境 的 胁迫 ,此 
外 ,分 泌 的 胞 外 水 解 酶 能 水 解 动物 和 植物 的 残骸 和 
土壤 中 一 些 复 杂 的 有 机 化 合 物 ”” ,这 可 能 是 部 善 
库 木 塔 格 沙漠 和 塔克拉玛干 沙漠 等 这 类 流动 沙漠 
中 优势 菌 门 是 放 线 菌 门 的 原因 。 

LEfSe 分 析 表 明 ,在 门 水 平 上 Campilobacterota 
(如 Collinsella) 和 NB1-j 在 表层 样本 中 显著 富 
(P<0.05) , Sumerlaeota 在 深层 样本 中 显著 富 集 。 
Campilobacterota 是 动物 肠 道 的 优势 菌 门 ,并 且 发 现 
高 温 环境 会 使 Campilobacterota FY fal FERS IN. 5 
外 ,Collinsella 被 发 现存 活 于 人 类 上 肠 道路 。 因 此 ， 
推测 以 Collinsella 为 代表 的 Campilobacterota 在 表层 
样本 中 出 现 富 集 的 原因 可 能 是 由 于 采样 点 位 于 部 
善 库 木 塔 格 沙漠 北 缘 , 靠 近 人 类 居住 地 ,在 风 的 搬 
运 下 动物 或 人 的 凑 便 被 传播 于 沙漠 的 表层 ,从 而 采 
样 时 被 检测 出 来 。Sumerlaeota 第 一 次 被 报道 是 在 
被 水 淹 的 缺 氧 水 稻 土 壤 中 ” 。 另 外 ,Sumerlaeota 在 


EE ARK TER i ROY TEES CC 
和 人 工 湿 地 ”等 各 种 极端 环境 中 也 被 发 现 。 研 究 
表明 , Sumerlaeota 能够 降解 难 降解 的 有 机 物 作 为 碳 
源 ,还 具有 硝酸 盐 和 硫酸 盐 的 异化 还 原 .固氮 ERE 
和 有 机 磷 矿 化 的 能 力 ,使 其 在 各 种 恶劣 条 件 下 都 能 
苗 壮 成 长 ,此 外 研究 还 揭示 了 Sumerlaeota 的 成 员 是 
从 兼 性 厌 氧 的 共同 祖先 进化 而 来 的 拉 。 以 上 研究 
表明 , Sumerlaeota 能够 在 缺 氧 或 者 氧气 稀少 的 环境 
中 生存 。 另 外 ,部 善 库 木 塔 格 沙漠 北 缘 深 层 样本 中 
发 现 Noviherbaspirillum 宣 集 显著 ,推测 Noviherbaspi- 
rillum 能够 在 缺 氧 或 天 氧 环境 中 生存 。 造 成 都 善 库 
木 塔 格 沙 漠 北 绿 表层 和 深层 样本 的 细菌 群落 富 
种 类 在 门 水 平 上 存在 差异 ,原因 一 方面 可 能 是 不 同 
深度 的 含 氧 量 不 同 , 另 一 方面 可 能 是 不 同 深度 土壤 
提供 给 微生物 的 营养 物 不 同 。 已 有 研究 表明 , 阿 卡 
塔 玛 沙漠 不 同 深度 的 沙子 蕴含 着 独特 的 微生物 菌 
门 ,这 表明 不 同 深度 的 土壤 中 蕴含 的 微生物 具有 特 
定 的 偏好 性 ” 。 


4 # it 


(1) 部 善 库 木 塔 格 沙 漳 北 缘 沙 土 样本 呈 碱 性 。 
表层 的 不 同样 本 pH、EC、TC TOCH WC 35] 5 mH 
异 ; 深 层 不 同样 本 的 pH、TC、TN、TOC 和 WC 均 呈 显 
车 差异 。 

(2) 部 善 库 木 塔 格 沙漠 北 缘 优 势 菌 门 是 放 线 菌 
门 (Actinobacteria) 变形 菌 门 Proteobacteria) , EE fal 
门 (Firmicutes); 优势 属 是 unclassified_f_Micrococca- 
ceae „Bacillus ,Sphingomonas fl Escherichia-Shigella « 

(3) 部 善 库 木 塔 格 沙漠 北 缘 不 同样 点 之 间 细 菌 
群落 存在 显 车 差异 ;表层 和 深层 样本 之 间 的 细菌 多 
样 性 差异 不 显著 。 

(4) 属 水 平 上 表层 样本 有 17 个 显著 富 集 的 物 
种 ,深层 样本 有 14 个 显著 富 集 的 物种 。 

(5) WC、TC、TN 和 pH 均 显 著 影 响 痢 善 库 木 塔 
格 沙漠 北 缘 表层 细菌 群落 结构 ,TOC 显著 影响 深层 
样本 的 细菌 群落 结构 。 
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Study on bacterial community structure and influencing factors in the northern 
margin of the Shanshan Kumtag Desert 


LI Juan, LIU Yang", LIU Guangxiu™, CHENG Liang', 
GUO Qingyun, ZHANG Wei", ZHANG Gaosen 


(1. Institute of Plant Protection, Qinghai Academy of Agriculture and Forestry Science, Qinghai University, Xining 
810016, Qinghai, China; 2. Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, 
Lanzhou 730000, Gansu, China; 3. Key Laboratory of Extreme Environmental Microbial Resources and Engineering, 
Lanzhou 730000, Gansu, China) 


Abstract: The Shanshan Kumtag Desert, which is connected to the city, contains rich and unique microbial 
resources. To explore the microbial community structure and its relationship with environmental factors in the 
Shanshan Kumtag Desert, surface (0-5 cm) and deep (45-50 cm) layer sand samples were collected from four 
points along the northern margin of the desert. Subsequently, the collected samples were analyzed to determine 
the water content (WC), pH, total carbon (TC), total nitrogen (TN), total organic carbon (TOC), and electrical 
conductivity. High-throughput sequencing technology was used to analyze the bacterial community structure of 
the sand samples, and Spearman correlation was used to analyze the correlation between environmental factors 
and the bacterial community structure. The results showed that the collected sand samples were alkaline, and the 
dominant bacterial phyla were Actinobacteria, Firmicutes, Proteobacteria, and Bacteroidetes. Among these, the 
dominant genera were unclassified f Micrococcaceae, Bacillus, Sphingomonas, Escherichia- Shigella, and 
Microbacterium. Significant differences in the bacterial community structure were observed among the different 
points along the northern margin of the Shanshan Kumtag Desert (P « 0.05). However, no significant difference 
in species diversity was observed between the surface and deep layer samples. Furthermore, WC, TC, TN, and 
pH significantly affected the bacterial community structure in the surface samples (P « 0.05), while TOC 
significantly affected the bacterial community structure in the deep samples (P « 0.05). This study provides a 
theoretical basis for exploring microbial resources in the Shanshan Kumtag Desert and similar habitats. 
Keywords: Shanshan Kumtag Desert; high-throughput sequencing; bacterial community structure; environmen- 
tal factors; physicochemical property 


